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Abstract 

The vanishing of the divergence of the total stress tensor (magnetic plus kinetic) in a neighborhood of an equilib- 
rium plasma containing a toroidal surface of discontinuity gives boundary and jump conditions that strongly constrain 
allowable continuations of the magnetic field across the surface. The boundary conditions allow the magnetic fields on 
either side of the discontinuity surface to be described by surface magnetic potentials, reducing the continuation prob- 
lem to that of solving a Hamilton-Jacobi equation. The characteristics of this equation obey Hamiltonian equations of 
motion, and a necessary condition for the existence of a continued field across a general toroidal surface is that there 
exist invariant tori in the phase space of this Hamiltonian system. It is argued from the Birkhofl" theorem that existence 
of such an invariant torus is also, in general, sufficient for continuation to be possible. An important corollary is that 
the rotational transform of the continued field on a surface of discontinuity must, generically, be irrational. 
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1. Introduction 

This paper is concerned with plasma confinement in 
nonaxisymmetric toroidal magnetic fields, such as oc- 
cur in fusion experiments. As is typical in plasma equi- 
librium theory we assume the plasma to be static (i.e. 
mass flow is negligible) and the ion and electron Lar- 
mor radii to be negligible. There is no minimum scale 
length in this model and discontinuities are in principle 
possible, provided the net force density at each point 
r in the plasma vanishes (the force balance condition). 
Also, a corollary of the assumed flowless state is that 
electric fields are also negligible. 

We assume the kinetic stress is described by the 
isotropic pressure tensor P(r)l, where I is the unit 
dyadic. Adding the electromagnetic stress, the total 
stress tensor is 

T = PI+— l-BB, (1) 

where B = Bsi(r)/ yfJiQ, Bsi being the magnetic field in 
SI units and /io the permeability of free space. The force 
balance condition is then 1 1] 

V-T = 0. (2) 

If P and B are differentiable, Eq. (j2]) may be written in 
the more usual form Vf - (VxB)xB; but at surfaces of 
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discontinuity P and B are step functions and V X B con- 
tains a Dirac delta function component, corresponding 
to a sheet current at the discontinuity. Then (V X B) X B 
is not well defined and it is better to regard Eq. in- 
terpreted in the weak sense of distribution theory, as the 
fundamental form of the force balance relation. 

Magnetic field lines can be regarded as trajectories 
of a 1 5 degree-of-freedom Hamiltonian dynamical sys- 
tem [2, e.g.], with the toroidal angle ^ taken to be the 
analog of time and the half degree of freedom implying 
the system is not, in general, autonomous — the Hamil- 
tonian depends explicitly on ^. An important excep- 
tional case is that of axisymmetry, as in an ideal toka- 
mak, in which case the field-line Hamiltonian is inte- 
grable and all field lines lie on invariant tori foliating 
the plasma volume. It is then consistent with Eq. (|2]) to 
assume P and B to be arbitrarily smooth functions, with 
the isosurfaces of P coinciding with the invariant tori. 

In the theory of magnetically confinement plasmas 
these invariant tori are called magnetic surfaces, and 
each is characterized by a field-line winding number 
called the rotational transform and denoted t, being the 
asymptotic limit of the number of poloidal transits di- 
vided by the number of toroidal transits as the length of 
a field-line segment is taken to infinity. 

However, in stellarators and real tokamaks axisym- 
metry is (to a greater or lesser extent) broken and the 
magnetic field is nonintegrable. Then, while magnetic 
surfaces of sufficiently irrational *, may exist [by the 
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Kolmogorov-Arnold-Moser (KAM) theorem for small 
departures from axisymmetry, or by good design in sys- 
tems far from axisymmetry], the intervening magnetic 
surfaces are broken and replaced by a fractal structure of 
island chains, which include periodic orbits and chaotic 
regions. 

In this case the problem of finding P and B glob- 
ally consistent with Eq. (|2]) is notoriously difficult lO, 
and it is an open question as to whether solutions with 
globally continuous but nonconstant P exist. (In the lin- 
ear force-free case of constant P throughout the plasma, 
when B obeys the Beltrami equation, V X B = //B, with 
ji - const, solutions are known to exist under appro- 
priate boundary conditions 1^,31, irrespective of the in- 
tegrability of B.) However, by appeal to KAM theory 
it has been shown ||6l that discontinuous solutions with 
piecewise constant P exist for small departures from ax- 
isymmetry, and recent work [Z, i8J holds out promise 
that variational methods can be used to construct numer- 
ical solutions with a stepped pressure profile for systems 
far from axisymmetry. 

In the present paper we consider a simpler, more local 
problem: Consider a prescribed torus iS and assume P 
and B differentiable in open regions bounded by S, with 
P constant on each side of S but discontinuous across 
it. Then, given the field B(r) on one side of S and i 
on the other side, find the field B(r) on the other side 
that is consistent with force balance in the weak sense 
discussed above. 

In Sec. [2] we show that the surface magnetic fields 
must be derivable from potential functions defined on 
either side of the surface, so that the prescribed B can be 
obtained by specifying a scalar function. The unknown 
surface potential on the other side can be determined by 
solving a Hamilton-Jacobi equation involving a Hamil- 
tonian, the pressure jump Hamiltonian. The configura- 
tion space on which this Hamiltonian is defined corre- 
sponds to a single flux surface, while the "momenta" 
correspond to possible magnetic fields on the same sur- 
face. Thus the pressure jump Hamiltonian is quite dis- 
tinct from the magnetic field line Hamiltonian described 
above, which is defined globally. (These Hamiltonians 
are compared and contrasted in Appendix B. The only 
relevance of the magnetic field line Hamiltonian to the 
problem posed in this paper is that S must be an invari- 
ant torus of the field line Hamiltonian dynamics.) 

Earlier, Berk et al. |9| introduced a Hamilton-Jacobi 
approach to investigate an individual flux surface's abil- 
ity to withstand pressure, but their scope was limited to 
the case of vanishing current field on one side of the sur- 
face and vanishing magnetic field on the other Kaiser 
and Salat IfTOl attempted to prove the existence of flux 



surfaces using only geometrical arguments (geodesic 
flow), but this approach assumed the same field simpli- 
fications as Ref. [9 1. By contrast our formulation allows 
nonzero current and magnetic field on both sides of S. 

The geodesic approach of Ref. 1 10] was introduced 
because of a worry that points on an invariant torus em- 
bedded in the phase space of the pressure jump Hamilto- 
nian may not map one-to-one to the configuration space 
of the torus, violating the condition that field lines can- 
not cross. Kaiser and Salat were also doubtful of the ap- 
plicability of the KAM theory (which had been alluded 
to in Ref. |i9J) to ensure the existence of some flux sur- 
faces under small deformations. 

In Sec. [3] we demonstrate the applicability of the 
Hamiltonian formulation by applying the Birkhoff theo- 
rem to prove the existence of a mapping that takes char- 
acteristics on invariant tori in the phase space of the 
pressure jump Hamiltonian to the magnetic field lines 
in configuration space, such that the topology of those 
characteristics is preserved. We then briefly discuss 
conditions and methods for finding such invariant tori, 
arguing that, while KAM theory is applicable to sys- 
tems close to axisymmetry, this is too restrictive to be of 
much practical use, whereas Greene's residue method is 
practical if not completely rigorous. 

One might be tempted to infer from this, and the 
fact that invariant tori in nonintegrable Hamiltonian sys- 
tems generically have irrational winding numbers, that 
( must be irrational on both sides of S. However, in the 
plasma equilibrium problem, the field line dynamics is 
not known a priori but rather to be found by solving 
for B self-consistently, with a boundary condition being 
that B be tangential on the two sides of S. Thus, it could 
be said that S is an invariant torus of the field line dy- 
namics by construction. It is one of the main aims of 
the present paper to show that force balance across S 
provides an independent, locally self-consistent reason 
for believing that t must be irrational on both sides of S 
in almost all nonaxisymmetric equilibria. 

2. Formulation 

2.1. Physical Derivation 

We are interested in weak solutions of the coupled 
pair of equations, Eq. (|2]i and the divergence-free mag- 
netic field condition. 



V-B = 



(3) 



in the neighborhood of a given surface of discontinuity. 

Much can be inferred from Eq. (|2]l; the details are in 
Appendix Appendix A but here we report the important 
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Figure 1 : (color online) A toroidal cut showing the cross section of 
the plasma and the various definitions. The cross section is a circle 
for simplicity, but the theory will work for a surface S of any shape. 



results. When there is a discontinuity in the pressure 
across a surface S, and when the magnetic field is force 
free in the neighborhood of the surface, the following 
must apply on S, 



P - const . 
B-n = , 
(V X B) • n = , 



1 7 

-B- + P 
2 



= 



(4a) 
(4b) 
(4c) 

(4d) 



In Equations n is the normal to the surface, and 
the notation IjcJ refers to the difference between x on 
one side of the surface and x on the other side. Equa- 
tions (|4a]-|4c]) must hold on both sides of the surface 
[with a different constant on either side in the case of 
Eq. (4ai]. We refer to Equations (Hh collectively as the 



pressure jump conditions. Note the tangential compo- 
nent of the magnetic field is discontinuous. 

To describe toroidal fusion geometries we now con- 
sider the surface S to have a toroidal topology and use 
the coordinate system {0, Here, 9 and f are angle- 
like coordinates in the poloidal and toroidal dkections 
respectively. 

The field on the inner side (side closest to the mag- 
netic axis) of S is given by B , and the field on the outer 
side is given by B^, as shown in Fig.[T] Equation (4di 
then becomes 



B+2 _ ^ 2(P_ - P+) . 



(5) 



When Eq. (jSj) is written in terms of the covariant com- 
ponents of the magnetic fields, one finds that 

2AP = Yu [^"^1 - ^i^l] ' (6) 



where AP = P - P+ and g'^ = g'^{0,O are 
the contravariant metric coefficients defined on the 2- 
dimensional Riemannian manifold S, related to the co- 
variant coefficients |j6l by inverting the 2x2 matrix [gij], 



8 -8((/^' g --ged^^ 8' 



^8: 



(7) 



where ^[g - 8ee8ii ~ 8b( '■^^ determinant. We shall 
also need the 2-dimensional contravariant components 
of the surface magnetic field, defined by 



a -g «fl+^ 

= g'^^Be + g((B^ . 



(8a) 
(8b) 



The torus S is embedded in an ambient, 3- 
dimensional Euclidean space with its own metric, G'-' 
say. If {6,(,iff) are curvilinear coordinates in this ambi- 
ent space such that tj/ - const on S, then the covariant 
components G,j, j 6 {6, are identical to gij on S, 
but the contravariant components are different. 



(9) 



2.2. Hamiltonian Treatment 

The keys to the Hamiltonian treatment are Equa- 



tions (4b I and (4c i, together they imply that on S, 

dgB^-d(B^=0. (10) 



Eq. (10 1 is implicitly satisfied if the field components 
are written as 



(11) 



where the two scalar functions /"^(O, ^) - JB"^ • dl are 
referred to as surface potentials. The surface potentials 
and define fully the fields on the inner and outer 
sides of S (B^ and B^) respectively 

Now with Equations ( [4b| and ( pcj i being implicitly 
satisfied, the pressure jump conditions reduce to the sin- 



gle condition, Eq. ( 4d i 



To this end we substitute Eq. (Ill into Eq. m« to give 



2AP^ 8''[derda^-dgf-dif-] . (12) 

The goal is, given a surface (which defines g'^) and a 
field on one side of the surface (which defines, say / ) 
the goal is to find the potential on the other side of the 
surface (/^) (See Fig. |2|. If can be found that has 
continuous second partial derivatives [/^ e C^(P^)], the 
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Figure 2: (color online) A depiction of a toroidal segment of the sur- 
face S to demonstrate the properties required for the Hamiltonian. 



magnetic field will satisfy force balance and lie on the 
surface, satisfying Equations Q. 

The problem is symmetric, we can either, given 
find (work inside-out), or given find (work 
outside-in). So we refer to the unknown surface field 
potential simply as /. Then the pressure jump condition 
becomes a problem of calculating / from the equation 



canonical pairs (6, pg) and (f , ) through 

Pe = def ^ Bg , = d^f = . (16) 

If / exists, then the Hamiltonian orbit lies on an in- 
variant torus in phase space. Such an orbit we refer to 
as regular. 

Treating the problem as Hamiltonian, one is able to 
utilize tools that have been developed for determina- 
tion of integrability in Hamiltonian systems to investi- 
gate whether a solution (/) can be found that satisfies 
force balance. The intent of this paper is to provide a 
consistent explanation to prove that inferences from the 
Hamiltonian system are acceptable. 

Determination of the existence of a surface potential 
is sufficient in the sense that it dictates definite existence 
criteria for the data supplied-surface shape, rotational 
transform, prescribed magnetic field on one side and 
pressure. However, on the scale of the entire plasma 
system it is only a necessary condition in that many 
other sources of destruction may be present throughout 
the plasma. These perturbations are generated outside 
the domain of the pressure jump Hamiltonian as it is 
only defined within an infinitesimal region of the sur- 
face. 



H(0,(,dgf,dj) ^ AP 



(13) 



Eq. ( 13 I is a partial differential equation for /. More 



specifically, Eq. (13 1 is a time independent Hamilton- 
Jacobi Equation with 5,/ = pi and Hamiltonian 



H(0,(,pg,P()^g'JpiPj + V(6,O^ 



(14) 



where there is an implicit sum of i,j e {6,(}, and 
V(0,^) = g'^dif^djf^ where has been arbitrar- 
ily chosen as the given potential. The Hamiltonian in 
Eq. ( 14 1 we refer to as the pressure jump Hamiltonian. 



The solution to this autonomous Hamiltonian system 
can be found by solving the corresponding characteris- 
tic equations 11 li 



Pe 



dH_ 
dpg ' 
_dH 

30 ' 
dH 

9pc ' 
_dH_ 



(15a) 
(15b) 
(15c) 

(15d) 



As a Hamiltonian system, AP is identified with the 
energy and / is a type two generating function (Hamil- 
ton's characteristic function) considered to generate the 



2.3. Reduction to a DOF system 

To simplify computation of Hamiltonian orbits, we 
condense these equations by dividing Equations ( fTS) by 
Eq. ( 15c I to trivialize the third of Hamilton's equations 
and make the toroidal angle-like coordinate the "time" 
variable. The division requires that ^ 5^ 0, the implica- 
tions of this being addressed in Sec. 



3.1 



The reduced equations provide more physically rele- 
vant equations of motion. The first describes the path of 
the Hamiltonian trajectory through configuration space: 



do _ g^^pe + g'^^P( _ If 
d^~ g'^Pe + g^^Pi'B^ 



(17) 



where Equations (|8]l have been utilized. Eq. ( 17 1 is the 
equation of a field line. Any solutions to this charac- 
teristic equation of the pressure jump Hamiltonian may 
correspond to field lines of the magnetic field on S. 
The second equation of motion is 



dpi _ djg'JpiPj- djV 
d( ~ gO^pg + giip^ 



(pi = B/) 



(18) 



for / € {0,^}. Eq. (18 1 shows that the canonical mo- 
mentum gives the covariant components of the magnetic 
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Physical Quantities and their Hamihonian Equivalents 



AP 


Pressure jump 


E 


Energy 


f 


Surface potential 


S 


Action / Hamilton's principal function 




Curvilinear coordinates 


q 


Generalized coordinates 




Covariant components of magnetic field 


p 


Generalized momenta 


(©,<!)) 


Straight field line coordinates 


Q 


(Action) Angle coordinates 


t 


Rotational transform 


w 


Winding number (angular frequency) 



Table 1 : A table summarizing the physical interpretations of Hamiltonian quantities in the problem. 



field along a Hamiltonian trajectory for solutions with 
the required rotational transform. 
The final equation of motion, Eq. 



15di, can be solved 



implicitly by the first two using the fact that the en- 
ergy of this Hamiltonian system is conserved, i.e. the 
pressure jump is constant along the flux surface because 
neither region adjacent to our flux surface interface can 
support a pressure gradient. This means can be writ- 
ten as a function = p^{9,pg,(;AP). When this is 
substituted into Equations (|15a[l and (15bi, the entire 



system can be solved within two differential equations. 
However, inversion of the Hamiltonian brings about an 
arbitrary sign, which for simplicity, we choose to be 
positive, discussing the choice further in Sec. |3.1[ 

The system is now condensed into the differential 
system 



d0 
dpe 



u{e,(,pB\AP) 

v{0,^,pe-AP) 



(19a) 
(19b) 



One method of identifying the existence of / is to calcu- 
late the Hamiltonian trajectories. If a solution to this 1 1 
degree of freedom system can be found that lies on an 
invariant torus in (6,pg,() space with the correct wind- 
ing number, the trajectory, when projected onto the 3D 
geometric torus, coincides with the field lines that lie on 
that surface. 

The solution of the pressure jump Hamiltonian is in 
general not unique. It may be that many Hamiltonian 
trajectories are regular, and so are physical candidates 
for a field to satisfy force balance. To make the solution 
unique one can require the corresponding field line have 
a certain rotational transform, defined as 



lim 



A9 
A? 



(20) 



The pressure jump conditions. Equations Q, allow a 
jump in i across S, although there is some evidence 
to suggest that a jump in rotational transform is unsta- 
ble d. 



Given the non-intersection of phase-space character- 
istics for a well-defined Hamiltonian system |131 and 
the fact that the phase-space characteristics are confined 
to a topological torus, it is understood that the rotational 
transform embodies a topological invariant of the phase- 



space characteristic itself. In Sec. 3.2 the existence of a 
map between phase-space characteristics and the mag- 
netic field lines in configuration space that preserves 
characteristic topology will be shown. It is clear that 
since this map preserves characteristic topology, it also 
preserves the rotational transform. 

2.4. Action Angle Coordinates 

In action angle coordinates the canonical momenta 
are constant, so that the equations of motion are trivial. 
The equations of motion are, after an appropriate choice 
of initial conditions. 



— WQt , 



where wq and wo are constant. Thus, 



(21) 



(22) 



This coordinate system is, in fusion research, known as 
straight field line coordinates, as the magnetic field ap- 
pears as a straight line in these coordinates. Such co- 
ordinates are helpful because the rotational transform, 
usually a quantity that requires integration along the en- 
tire length of a field line (often infinitely long) is now 
explicit in the equations of motion: 



AO d@ 
lim — = — 
Af^oo A( d(i> 



Wo 



(23) 



Such coordinates can be found when / (which is 
the generating function, cf. Eq. ( [T6| ) is separable 
in the configuration coordinates, in this case the cor- 
responding Hamiltonian trajectory lies on an invariant 
torus |14|. 

The most general form for / given its definition in 
Eq. ( [TT] ) is 



(24) 



6 



where / and G are constants, and f{6, is a function 
periodic in and ^. The transformation to straight field 
line coordinates can be accomplished via the transfor- 
mation 



= 0- 
In the new coordinates. 



mo 



(25a) 
(25b) 



/ = /© + GO , (26) 
and thus the magnetic field components are constant. 

3. Analytical Concerns 

3.1. Ambiguity of Sign 

The pressure jump Hamiltonian is a constraint on the 
square of the magnetic field, so it is expected that there 
are two magnetic fields that would satisfy force balance. 
This arbitrariness is made explicit when one inverts the 
Hamiltonian to find = Pi:{6, pe, ^; AP) in an effort to 
reduce the phase space. When completing the square 
one has the expression 



Pe + g^^Pi: = ±A , 



(27) 



where A = (/^ [lAP + V(6, - ^>fl)+/^ pI^'^- The 
left hand side of Eq. ( 27 1 is the result of Eq. ( |15c[ l, i.e. 



dt 



= = ±A , 



(28) 



showing that the two solutions correspond to magnetic 
fields with opposite toroidal direction. 

The information on the toroidal direction of the given 
field is similarly lost within the pressure jump Hamilto- 
nian so the choice of sign must be made to reflect the 
initial conditions. It is expected that physical configu- 
rations would require a field to be in the same toroidal 
direction on either side of an infinitely thin flux surface, 
otherwise there would be a very strong current sheet, 
susceptible to a tearing instability 1 15 1. 

The choice of sign has also been shown in Bruno and 
Laurence to be equivalent to choosing the sign of the 
rate of change of toroidal flux |6|. 

When reducing the phase space of the Hamiltonian 
system, Hamilton's equations were divided by ^, thus 
the condition in which ^ = was necessarily lost. A 
field line that does not extend toroidally corresponds to 
a field configuration of infinite rotational transform-a 
situation we ignore as it will not ergodically cover the 
surface, and thus wiU never act as a flux surface. 



3.2. Birkhoff Theorem 

Kaiser and Salat 1 10] solve the pressure jump discon- 
tinuity problem purely in configuration space, that is, a 
solution to force balance on the surface is sought that 
corresponds directly to geodesies covering a 3D torus. 
Such an approach is limited to situations where the field 
within the plasma volume is zero. However, Kaiser and 
Salat felt obliged to use this geodesic method because of 
concerns regarding the physical significance of Hamil- 
tonian trajectories in phase space. 

Kaiser and Salat's gravamen against the Hamiltonian 
formulation can be stated as the following: Suppose 
a solution to the pressure jump Hamiltonian system is 
found. This will correspond to a Hamiltonian orbit that 
lies in a four dimensional phase space. The actual field 
line however exists on the two dimensional torus em- 
bedded in Euclidean 3-space, i.e. configuration space. 
We must project a four dimensional phase space trajec- 
tory to the two dimensional configuration space-is it not 
possible that the projected trajectory intersects itself? 

Assuming the magnetic field is nowhere zero, such 
intersections would make it impossible to interpret the 
projection as a physical field line. 

However, we will now prove that such crossings can- 
not occur, via a direct application of the Birkhoff theo- 



rem to Equations ( 19 1. It will thus be demonstrated that 
the existence intrinsic to the Hamiltonian formulation 
is sufficient to imply that the corresponding field line is 
consistent with Equations Q. 

Our system is a 1 ^ degree of freedom Hamiltonian 
whose trajectories define a 2D area preserving map by 
integrating Equations (19 1. The mappings of interest. 



those generated by a trajectory that Ues on an invariant 
surface, are also twist maps as the metric is positive def- 
inite [dtl{dpdp.H) > 0] |16|. 

Consider the phase space variables {q,p) in a 2 di- 
mensional area preserving twist map, the Birkhoff The- 
orem states that, for a rotational invariant circle, Iil7i 



P = Y(q) , (29) 
where y is a Lipschitz function on i.e. Y satisfies 
\Y(x)-Y(y)\ 



sup 



\x-y\ 



< c 



(30) 



for some bounded constant C. In this case we can write 
the phase space mapping generated by the Hamiltonian 



(q',p')^T(q,Y(q)) 



(31) 
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Let us consider the operator n that is the projection of 
the phase space trajectory onto configuration space. 



then 



7i{q, p) = q. 



q' = n{T[q, Y(q)]) = a(q) . 



(32) 



(33) 



Thus, as r is a homeomorphism and Y Lipschitz, a is 
also a homeomorphism |,.17J . The injective nature of a 
homeomorphism imphes there will be no crossings un- 
der the mapping n. 

Strictly, this is only true for a two dimensional sys- 
tem because the Birkhoff theorem apphes only for a 2D 
phase space. Some limited higher dimensional results 
have been found [18J . 



This means a homeomorphic mapping like Eq. ( 29 1 



can be generated to define completely the evolution of 
the system, and the above proof applies. Thus, when 
mapping the Hamiltonian trajectories to the 2D torus in 
configuration space no crossings are possible. 



whether the field line lies on an invariant surface. The 
most commonly used tool to accomplish this is Green's 
residue criterion |19|. Despite only partial mathemat- 
ical justification, 1 20 1 has proved helpful in investiga- 
tions of this kind 12111221 12I1. 

A more computationally expensive but perhaps more 
physically relevant method would be the calculation of 
the analyticity width Il24l [20 1. Technically, it repre- 
sents the domain of analyticity of the transformation to 



action-angle coordinates in Equations ( 25 1 



There has been some progress in non-Hamiltonian 
approaches to investigate existence in problems like the 
pressure jump Hamiltonian. Kaiser and Salat computa- 
tionally found evidence for KAM-like behavior in their 
purely configurational treatment. They developed their 
own purely configurational theorem to determine the ex- 
tent of deformation such that no two field lines can lie 
on the torus and not intersect, which they termed the 
"big bump" criterion IfTOl . 

4. Conclusion 



3.3. Existence of Invariant Tori 

With the acceptance that solutions garnered from the 
pressure jump Hamiltonian map homeomorphically to 
field lines on the 3D toroidal surface, we can use the 
existence of the surface potential to declare that the cor- 
responding field line satisfies Equations Q. 

When a Hamiltonian trajectory can be transformed to 
action angle coordinates, the Hamiltonian orbit lies on 
an invariant torus; and after the mapping to configura- 
tion space the field line will lie on the surface and not 
intersect itself. Conversely, if the conditions are such 
that the Hamiltonian trajectory is chaotic, it does not lie 
on an invariant torus; after the mapping the field line 
wiU not lie on the given flux surface (B • n = is not 
satisfied) and no physically consistent field can exist. 

Of great interest then, is the knowledge of whether an 
invariant torus of the Hamiltonian exists or not. There 
are various tools to use. Previous treatments of the prob- 
lem have suggested using the KAM theory to prove the 
persistence of flux surfaces under small perturbation in 
the Hamiltonian |9 |. In the pressure jump Hamiltonian, 
it is not clear in advance whether a given trajectory is an 
invariant surface, and so the question of its persistence 
under perturbations is not helpful. 

Practically, especially in fusion devices, the surface 
would be strongly perturbed and it would not be clear if 
the flux surface with a given rotational transform would 
exist. Thus existence must be determined a posteriori, 
that is, after the field line has been found, determine then 



By applying force balance to neighboring regions 
of plasma of finitely different pressure separated by a 
infinitesimally thin flux surface, three general criteria 
were found that must be satisfied in order for the sur- 
face to be a flux surface. The three criteria were com- 
bined into one, referred to as the pressure jump criterion 
by condition of the existence of a surface potential /, a 
scalar function defined on the surface that is related to 
the magnetic field. 

A Hamiltonian-Jacobi construction of the pressure 
jump criterion was introduced, and the surface potential 
was found to play the role of the generating function 
important in Hamiltonian dynamics. 

The main analytical problem that faced the treatment 
was solved, namely the sufficiency of an existence crite- 
rion on / to determine one-to-one that the correspond- 
ing field line satisfied the pressure jump criterion. This 
was resolved by appealing to the Birkhoff theorem when 
the Hamiltonian system is integrated to form a two di- 
mensional area preserving map. The application of the 
KAM theory to Hamiltonians of this type has been ques- 
tioned in other papers, so practical avenues of investiga- 
tion other than the KAM theorem were suggested. 
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Appendix A. Pressure Jump Condition 

Consider a surface S, defined as n = where n 
is the distance from S. For the purposes of deriv- 
ing the pressure jump condition, the definitions of B"^ 
are extended radially so they are both smooth every- 
where (i.e. not by themselves discontinuous) and re- 
main bounded in both volumes. The discontinuous pres- 
sure profile P(r) and the corresponding discontinuous 
complete magnetic field B(r) can then be expressed as 

P = P^{r)h(-n) + P+{r)hin) , (A.l) 
B = B_(r)/!(-«) + B+(r)/z(n) , (A.2) 

where h is the unit Heaviside step function, which 
causes the discontinuity in B, as both B'^ are themselves 
not discontinuous. Such a form for P and B gives 

VP±=j^xBS (A.3) 

where jj_ = V X Bj- are the associated currents. Dotting 



Eq. (A.3 1 with B-t gives B-t • WP± = 0, which implies 
the pressure is constant along a magnetic field line. As a 
flux surface is composed of a single magnetic field line, 
we have the condition that on both sides of the surface 



P± — const . 
Substitution of Eq. (|A.2|i into V ■ B = gives 



(A.4) 



V ■ B = (A.5) 
n- m S(n) + V ■ BJi(-n) + V ■ B+/i(n) , (A.6) 

where n is the unit normal to S and where Hjc]] is the 
jump of X across the interface, [jc]] = x+ - x^. The 
divergence can only be zero if 



IB J = n . IBl = , 



(A.7) 



i.e., the normal component of the magnetic field must 
be continuous. 

Similarly, when Equations ( |A.l - A.2 1 are substituted 
into the stress tensor in Eq. 



PI -B^-BB 



. 



(A.8) 



Substituting Eq. ( A.7 i into the above and dotting with n 
gives the condition 
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P+ 



B„ PJ = . 



(A.9) 



removing the middle equality gives the pressure jump 
condition. 



As B„ is continuous, Eq. ( A.9 1 can be written as 

|(n X B)^] = -2 IPl , (A. 10) 

so long as IPJ + 0, 



In X BJ 



Crossing Eq. ( A.8 1 with n gives 

B„ In X Bl = 



(A. 11) 



(A. 12) 



Combining Eq. (A. 11 1 with Eq. (A. 12 1 implies 

n-B±=0. (A.13) 

Thus, if there is a pressure discontinuity across a sur- 
face, the field lines must lie on that surface. 

The plasma in the neighborhoods either side of the 
surface is assumed to be force free, so that VP = 0. 



Eq. (A.3 I then implies that jj_ is parallel to B-t, which 



on comparison with Eq. (A.13 i implies 
n.j^ =n.(VxB±) = 



(A. 14) 



Thus, if there is a pressure discontinuity across a sur- 
face, the curl of the magnetic field must be parallel to 
the surface at all points on the surface. 

Appendix B. Comparison of Hamiltonians 



In this section we compare the magnetic field line 
Hamiltonian system often used in the literature to the 
pressure jump Hamiltonian system introduced in this 
paper. 

The magnetic field in the plasma volume be written 

asina 



B = V(/r X Ve -H X V;^ 



(B.l) 



Using the equation of the field line dridt - B(r), one 
finds la 



d( dt// 
dip dx 

~d^ ~ ~de 



(B.l) 
(B.3) 
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Magnetic field line Hamiltonian 



Pressure jump Hamiltonian 



P implicit 

Defined throughout plasma volume 

Phase space and time a representation of Euclidean 

3-space 

All orbits are field lines 



P explicit 

Defined on a given magnetic surface 

Phase space a combination of a Riemannian 2-space 

- magnetic surface (configuration space) and field 

components (momentum space) 

Orbits are not field lines, projections of regular orbits 

with specified * onto configuration space are field 

lines 



Table B.2: A table summarizing the differences between the magnetic field line Hamiltonian and the pressure jump Hamiltonian. 



which are of the form of Hamilton's equations. Thus 
each field line can be described as the solution of the 
magnetic field line Hamiltonian in ^ phase space 
with the poloidal angle and toroidal flux function 41 
as canonical variables. 

Trajectories of this Hamiltonian correspond directly 
to magnetic field lines within the plasma. Trajectories 
that lie on invariant tori correspond to field lines that 
draw out flux surfaces. Trajectories that are chaotic cor- 
respond to chaotic field lines the compose the chaotic 
regions of the plasma. 

The pressure enters implicitly into the magnetic field 
line Hamiltonian system, exciting currents that deter- 
mine the flux functions. The field line dynamics feeds 
back into the determination of the pressure as a function 
of position. 

In contrast, the pressure is explicit in the pressure 
jump Hamiltonian, but it only holds on a given toroidal 
surface. While the trajectories of this Hamiltonian sys- 
tem are not the field lines, the projection of the phase 
trajectories onto the 3D torus are field lines of the given 
rotational transform. If the surface potential can be 
found, then the surface is a flux surface. 

A summary of the differences of the Hamiltonian sys- 



tems in given in Table B.2 
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